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■ Abstract. Direct algebraic method of obtaining exact solutions to nonlinear PDE's is applied to certain set of 
nonlinear nonlocal evolutionary equations, including nonlinear telegraph equation, hyperbolic generalization 

, of Burgers equation and some spatially nonlocal hydrodynamic-type model. Special attention is paid to the 

■ construction of the kink-like and soliton-like solutions. 
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O ■ 1 Introduction 

^! In recent decades the problem of obtaining exact solutions of nonlinear evolutionary equa- 

(— I . 

^— > ■ tions has attracted attention of many experts in mathematical physics. The most fun- 

damental achievement in this area was development of the inverse scattering method pp. 
Unfortunately, this method is applicable to the relatively narrow class of the completely 

•i-H . 

^ ■ integrable equations, while for the majority of nonlinear evolutionary equations methods of 

I obtaining the general solution do not exist. Yet even the possibility of obtaining particular 

solutions to the nonlinear PDE's very often is considered in applications as a big success, 
because having the analytical solutions of a modelling system it is more easy to analyze it 
and to make its interpretation. Exact solutions are also widely used as a starting point for 
various asymptotic methods, for testing the numerical schemes and facilitating the stability 
analysis. 

For nonlinear PDE's, which are not integrable, exact solutions are usually obtained by 
means of the Group Theory Reduction |2j. Appreciating all advantages of the symmetry- 
based methods, we would like to pay attention to the fact that they are not fully universal, 
firstly, because they can be effectively used merely in case the PDE's possess some non-trivial 
symmetry, and, secondly, because it is very difficult within the frameworks of these methods 
to obtain any solution with the given properties. Note that much more efficient from this 
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point of view is the combination of the symmetry reduction with the methods of quahtative 
analysis. 

Aside from the symmetry reduction, there exist another group of methods enabhng to 
obtain exact solutions to nonlinear PDE's. They are based on choosing the proper trans- 
formation (or ansatz), simplifying the problem. As an example let us mention the famous 
Cole-Hopf transformation, which was originally used for the non-local linearization of Burg- 
ers equation. Later on employment of the similar transformations enabled Hirota to obtain 
the multi-soliton solutions of the completely integrable KdV equation without referring to 
the inverse scattering method [li . An intense development of the ansatz-based method in the 
following years was highly simulated by the fact that it proved to be effective for obtaining 
the particular exact solutions of the evolutionary PDE's that are not completely integrable. 
The essence of the ansatz-based method is well summarized in the recent work of E. Fan [Hj. 
Beside the imposing bibliography, it contains general formulation of the unified algebraic 
method, which will be presented in next section. Literally during last years the approach 
presented by E.Fan was used in the composition with the generalized Cole-Hopf anzatz, 
which made possible to obtain a series of exact solutions of nonlinear transport equation 

mini 

In this paper we present certain modification of the anzatz-based method presented in 
In accordance with the core of our interests, we use it to obtain exact solutions of 
hyperbolic modification of the non-linear transport equation, which takes into account the 
nonlocal effects. We concentrate on finding out the particular types of travelling wave (TW) 
solutions, namely, the soliton-like and kink-like solutions, but this approach can be easily 
used for searching out any exact solution which can be described as an algebraic combination 
of certain types of special functions. Next we consider the family of TW solutions for 
the spatially nonlocal hydrodynamic-type model. Using the qualitative analysis we state 
the existence of periodic and soliton-like TW solutions. Imposing some restrictions on the 
parameters we obtain the soliton-like solution in the analytical form. It turns out to be much 
more complicated than that obtained within the above mentioned anzatz-based methods. 
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2 Exact TW solutions to the hyperbolic modification 
of nonlinear transport equations 



2.1 Unified algebraic method and it's modifications 

The essence of the anzatz-based method or the direct unified algebraic method in E. Fan's 
terminology is based on the observation that the particular solutions of any system of 
PDE's 

H^iu', ul ui, ui^,...) i = l,2,...m, z/=l,2,..,n (1) 

which do not depend on {t, x) coordinates in explicit form, can be presented as a linear 
combination 

n 

«^ = E«>^(0 ^ = x + vt, (2) 

fi=0 

where unknown parameter, while function 0(0 - satisfies the equation 



0(O = ±. Ec.r- (3) 

Depending on the conditions posed on the parameters c^, solutions of the equation (0) are 
expressed by the elliptic Jacobi ( Weierstrass) , hyperbolic or trigonometric functions ^j. The 
properties of these functions are inherited by those solution of the initial system that can 
be presented in the form (j21). In fact this methodology is constructive and algorythmic if 
the functions H'^ arising in (^Q) are algebraic ones. One easily gets convinced that with 
this assumptions the substitution of into the initial system (P) gives the polynomial 
functions with respect to 4>^{^) and 0^0(0- Equating to zero the coefficients standing at 
the corresponding powers of these functions, we obtain nonlinear system of the algebraic 
equations which determines particular solutions of the initial system. 

The enforced version of the Fan's method was put forward recently in It was used 

for searching out the exact solutions of the equation 

ut + Auu^- Ku^^ = f{u). (4) 

More preciesly, there has been proposed the ansatz 

u= — i-zr , ^ = X + vt + xq, k = l,2,.., (5) 
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where z = J2 (^ti4''^{0j 0(0 is the function satisfying the Q- Owing to this combination, 
the multi-parameter famihes of the exact solutions were obtained in the situation when the 
pure Fan's methodology does not work j3]. 

On analyzing different versions of the anzatz-based method, one can conclude that effec- 
tiveness of their employment is based on the mere fact that the family of functions (f)^{^) 
and (j)'^<p{^) is closed with respect to the algebraic operations and differentiating. In view of 
this a quite natural generalization of the already mention ansatze would be as follows: 

where the function 0(^) still satisfies the equation Q, but, in contrast to (0), dependence 
between the functions / and g is not assumed from the very beginning. Effectiveness of the 
anzatz (jUI) is demonstrated in the following subsection. 



2.2 Exact travelling wave solution to the nonlinear hyperbolic 
equation 

Let us consider the following equation: 

TUtt + Auu^ + But - t^Uxx = f{u) = Y1 '^'^^"^ C^) 

where r. A, B, k are non-negative constants. For A = equation ((7j) coincides with the 
nonlinear telegraph equation; for A ^ it coincides with the hyperbolic generalization of 
Burgers equation, while for A = B = - with the nonlinear d'Alambert equation. The 
hyperbolic modifications of nonlinear transport equations arise in a natural way when the 
memory effects are taken into account |6 . 

The main goals of this paper are to obtain the exact soliton-like and kink-like solutions of 
equation (|7j) and to present the advantages of anzatz © compared to the already mentioned 
modification [4^ of the unified algebraic method. 

Assuming that the solitons and kinks can be expressed by powers of function sech{^), 
which is the particular solution of equation (jH)) and function sinh{C,) which appears in the 
odd derivatives of the function sech{^), we use the following ansatz: 

9{i) Tr'oC^sech^{ai) + sh{ai)Y.7=od.sech-{ai) 
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or, what is the same, 

Inserting anzatz (jHl) ( or Q) into (|7j) and executing of all necessary operations, we ob- 
tain an algebraic equation containing, respectively, functions sech^{a^), sech^ {a^)sh{aS,) 
or exp [/itt^]. Regarding them as the functionally independent ones and equating to zero 
corresponding coefficients, we go to the system of algebraic equations. We do not expose 
the details of these calculations since they are simple but cumbersome. To accomplish them 
we used the package of symbolic computation " Mathematica" . The results obtained are 
presented below. 

I. For arbitrary A, B and f{u) = Aq + \iu{t, x) + X2u{t, x)^ + X3u{t, x)^ function 

u[t,x] = - ; r-rrr (10) 

satisfies equation ((Tj) if the following conditions hold: 



Ao = ^^{BvA + hQ), Ai = j-^ [abohiBvQA + hQ^) + X^aoa.A 
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^2 = ^ [ablbUBvA + he) + AaA^e] , A = ^ [-2hablbl + \,A'] . 

Here and henceforth we use the notation h = a{v'^T — k), A = aibg — ao^i) © = 
ttibo + a^bi. Equation (fTUI) defines a kink-hke regime when 69^2 > and cto/^o 7^ a2/&2- 
Using the conditions we can express the unknown parameters from formula (fTUI) by 

the parameters characterizing equation ((7j), yet, in general case it is too cumbersome. It is 
much easy to do when a = 2y/^X^\^/v, oq = — ai = y^— A0/A2, bo = bi = 1. From these 
conditions we obtain the solution: 



u{t, x) = \J \2Tanh 



V— A0A2 , 
[x + vt) 



and conditions: v = A^(^^+V-^;^^_;-y3+i6.Ai .) ^ ^ For «: = 1, r = 0, A = 0, 5 = 1 

function ()10j) coincides with the solution obtained in j^. 

Another example of the kink-like solution defined by the formulae (fTUIl - pH) is as follows: 
2 

bo[l + exp{2 a 0\ 
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where 60 = (-A2 ± ^Ai - 4 Ai A3)/Ai, A = Xq = 0, a = -(2A2 + 3 60 Ai)/[45 &o], and 

±v/^(2A2 + 3&oAi) 

4Air + 46oA2(52 + 3Air) + 6gAi(252 + 9Air) ' 

II. For A = 0, arbitrary B and f{u) = Xo + Xi/2u{t, x)^ +Xiu{t, x)+X3/2u{t, x)^ +X2u{t, x)^ 
function 



u(t, x) 



(12) 



60 + 6ie'^(^+''*) 

satisfies the equation© providing that the following conditions hold: 
Ao = 2alalah/A'^, A1/2 = -2aoaia{3he + BvA)/A'^, 

Ai = 2a{h{3Q^ - A2) + BvAe)/A'^, A3/2 = -2bobia{5he + BvA)/A\ 

A2 = Qblblha/A\ 

This solution defines the solitary wave regime if 60^1 > 0, |ao|/|&o| = 1*^2 1/ 1^2 1 while 
ao/bo ^ 02/62- 

III. For B = 0, arbitrary A and f{u) = Xiu(t,x) + X3u(t,x) function 

~ -al-3 al a2 e< + 3 ai ai e^"? + 4 e^< ' ^ " ^ + (13) 
satisfies ((Tj), when Ai, A3 are positive and the parameters are as follows: 6 ai 02 = — ^AsTAi, 



a = a/Ai A3/A, V = ±y (A2/A3 + k) /r. This solution is always singular, because for arbi- 
trary values of the parameters the expression in the denominator of the formula ()13|) nullifies 
for some ^ G i?^. 

IV. Now let us consider the case A = B = 0. 
IV a. For f{u) = Aq + Xiu{t, x) + X2u(t, x)^ + X^uit, x)^ function 

^(^'^^ = bo + 26ie"(-+-*) + 6oe2"(-+-*) ' ^ 
satisfies equation ((Tj) when the following conditions hold: 

Ao = (^0 {2al bo — afbo — ao ai bi )a h/ A^ 

Al = {{albl + Aaoaibob^ + al{-Qbl + bl))ah/A'' 

(15) 

A2 = 3 6o (2ao&o "~ '^i^o^i — aobl)ah/A'^ 



A3 = -2 60 ibl-bl)ah/A 
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For ao 7^ 0, 60 7^ 0, |ai| + 7^ equation (fT^ defines the soliton-like solution. One of 
tlie parameters, contained in ()14j) can be cliosen arbitrarily, while the rest can be expressed, 
using the (fTHjl . as the functions of the parameters, defining equation (jZj). We omit doing 
this in the general case, but present one particular example. Thus, for 61 = 0, bo = 1 and 
arbitrary a we have the solution 

With ao = -A2/(3A3), ai = y^2(Al - Ai A3)/Ai, v = ±^[\, - Xl/{3 X3) + ^a^] /{r a^), 
with the additional condition Aq + Aiao + A2ao + Asajj = 0. 

VI b. For /(u) = Ai/2M(t, a;)-'^/^+AiM(t, x)+A3/2^t(i, a;)^/^+A2?^(t, x)^ we obtain the soliton-like 
solution 

""^^'"^^ " {boeM-+vt) + (26o + 46i)e"(-+''*) + bo)^ ^^^^ 

with 

Xy2 = -5ah/bi Xi = {12bo + 4:bi)ah/bi 

A3/2 = -(15 bl + 10 bo bi)a h/bi A2 = (6 fog 6^ + 6 bl)a h/bi. 

VI c. For f{u) = Ao + Ai/2M(t,x)^ + Xiu{t,x) + A3/2'u()f:, x)^ the localized wave pack 

(aoe2"(-+''*) + (2ao + 4ai)e"(-+-*)+ao)2 
«('^' ^) = (e^(x+..) + 1)4 (17) 

defines a solution of ((Tj) if the following conditions hold: 

Ao = 2ao(ao + ai) a h/ai A1/2 = (Sag + 6aoai) a h/ai 
Ai = (12ao + 4ai) a /i/ai A3/2 = — 5 a /i/ai. 

VI d. For f{u) = Xiu{t, x) + Xs/2u{t, x)t + A2'u(t, x)^ function 

" (aoeM-+-*)+2aie-(-+-*)+ao)2 ^^^^ 

with 

Ai=4a/i A3/2 = — 10 aia /i 

A2 = (6 a^ — 6 ao)Q; /i. 

defines the soliton-like solution of the equation ((7)). 
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VI e. Finally, let us consider equation 

TUtt - K,Uxx = Ao + XlU + X2U^ + Xsu^- (19) 

Inserting the anzatz u = 0(0; ^ = x + vt into the equation (|TI?|l . we obtain, after one 
integration, the following ODE: 



— = ± y Co + Ci M + C2 + C3 + C4 u^, (20) 

where Cq is an arbitrary constant, Ci = 2Xq/H, C2 = Xi/H, C3 = 2X2/(3 H), C4 = X^/{2H), 
H = Tv"^ — K. To this equation the classification given in (3j is applied: 

(a) if A2 = Ao = 0, then equation (fT^ possesses a soliton-like solution 

u{t, x) = ^J-2X,/Xssech{^X,/H0 Ai > 0, A3 < 0; (21) 

(b) if Ao = A2 = 0, then equation p9|) posesses a kink-like solution 

u{t, x) = ^J-Xi/Xs tanh ^^i/{2 H) ^ , Ai < 0, A3 > 0; (22) 

(c) if Ao = A3 = 0, then equation (fT^ possesses a soliton-like solution 

u{t, x) = — [3Ai/ (2A2)] sech 



Xi/Hi/2 



Ai > 0. (23) 



Presented above results enable us to state that the anzatze (jHI) and Q are effective and 
their employment gives the exact solutions in the situations when the ansatze suggested in 
jU 131 do not work. Note, that as a by-product we obtained a number of new exact solutions 
of non-linear transport and Burgers equations. These solutions can be easily extracted from 
the presented above formulae by simple substitution i? = l,r = (and also A = when 
it is necessary). Thus, the proposed modification of the unified algebraic method proves to 
be useful and its employment results in essential broadening the number of solutions of the 
given type (i.e. soliton-like and kink-like solutions), which can be obtained in analytic form. 
Yet, as it will be shown in the following section, none of the version of the anzatz-based 
method proposed by now is fully universal. 
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3 Periodic and soliton-like TW solutions of the nonlo- 
cal hydrodynamic-type model. 



In conclusion, let's analyze the family of TW solutions for the following system: 
ut + /Jp'^+Vx + <y [p'^^Px.. + 3(z/ + l)p>.p.. + u{u + l)p'^+V"" V'] = 0, 

, 2 n ^ ^ 

pt + p^Ux. = 0, 

where z/, /3, a are constants. The system arises in a natural when the balance equations 
for mass and momentum, taken in the hydrodynamic approximation, are closed by the 
dynamic equation of state, accounting for the short-ranged spatial non-locality [7j. In general 
case the answer on the existence of the periodic and soliton-like TW solutions is obtained 
by the methods of qualitative analysis, but under some additional conditions posed on the 
parameters we are able to present some of them in the analytic form, omitting the anzatz- 
based method (which doesn't work in this case). 

Let us consider the following family of invariant travelling wave solutions: 

u = U{uj), p = R{uj), uj = x-Dt. (25) 

Inserting the anzatz ()25j) into the second equation of system ()24j) . we obtain the first integral 
U = Ci-D/R and the system of ODE's 



dR _ Y 

duj 

d£ = {aR'+^y^ {ER - [D^ + /5i?'^+V(z/ + 2) + a{u + 1)R''+'Y^]} . 



(26) 



In accordance with the asymptotic conditions lim^^+oo f^(t^) = 0, lim^^^+oo -R(t^) = -Ri > 
0, f/(+oo) = 0, -R(+oo) = -Ri > 0, we assume henceforth that Ci = D/Ri and E = 
D^/R^+(3R'i+^/{u + 2). 

Dividing the second equation of system (|26|) by the first one and introducing new variable 
Z = Y^ = {d R/duj)"^, we get, after some algebraic manipulation, the linear inhomogeneous 
equation 

Z'{R) + 2[{u + 1) R\-'Z{R) = 2[eR-D^- PR'+yiu + 2)] /{a R'^^). (27) 

Solving this equation with respect to Z = Z{R) and next integrating the equation obtained 
after the substitution Z = [d R/duj)"^, we can express the solution of (j26j) as the following 
quadrature: 

f ±J^R^+''dR , , 

a; - cuo = / ^ 28) 



Unfortunately, the direct analysis of solution ()28p is very difficult. To realize what sort 
of solutions we deal with, the qualitative analysis of system is performed. It is evident, 
that all isolated critical points of system are located on the (horizontal) axis OR. They 
are determined by solutions of the algebraic equation P{R) = 13 R^^^ / [v + 2)—ER+D^ = 0. 
One of the roots of this equation coincides with Ri. Location of the second real root depends 
on relations between the parameters. If z/ > —2 and > PRi~^^, then there exists the second 
critical point A2{R2, 0) with R2 > Ri and the polynomial P{R) has the representation 

P{R) = {R- Ri){R- R2)'^iR). (29) 

For any u > —2 function \E'(-R) is positive whenever R > 0, since the function /?i?''+7(z/ + 2) 
is concave in this interval and has exactly two intersections with the line E R — D^. 

Analysis of system's (j26|) linearization matrix shows, that the critical points Ai{Ri, 0) is a 
saddle, while the critical point A2{R2, 0) is a center. Thus, system ()26|) has only such critical 
points, that are characteristic to the hamiltonian system. This circumstance suggests that 
there could exist a hamiltonian system equivalent to (j26|) In fact, the following statement 
holds. 

Lemma. // to introduce a new independent variable T, obeying the equation duo/dT = 
2(Ti?2('^+i), then system ^ can be written as a hamiltonian one 



dR _ 9 p2(i/+l) _ dH_ 

dT ^t^-n- — dv 



% = 2R^{ER- [Z}2 + 13 R-+^/[y + 2) + a{v + l)R^+^Y^]) = -||, 



(30) 



with 



By elementary checking one can get convinced that the function H is constant on phase 
trajectories of both systems (j26|) and (|30|). and since the integrating multiplier \E' = 2R'^, 
occurring in formula (jHUj) is positive for R > 0, then phase trajectories of systems and 
(jHUj) are similar in the right half-plane of the phase space [R, Y). Thus all the statements 
concerning the geometry of the phase trajectories of system (jHUj) lying in the right half-plane 
are applicable to corresponding solutions of system (j26j) . In particular, we immediately 
conclude that the critical point A2{R2, 0) remains a center when the nonlinear terms are 
added. This means that the initial system ()24|) possesses a one-parameter family of periodic 
solutions. If the right branches of the separatrices of the saddle Ai{Ri, 0) go to infinity (the 
stable branch when t —00 and the unstable branch when t +00), then the 
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domain of finite periodic motions is unlimited. Another possibility is connected with the 
existence of the homoclinic trajectories. In this case the initial system possesses localized 
soliton-like regimes. To answer the question on which of the above mentioned possibilities 
is realized in system (|26|). the behavior of the saddle separatices, lying to the right from the 
hne R = Ri, should be analyzed. We obtain the equation for saddle separatices by putting 
H = H{Ri, 0) = Hi in the LHS of the equation (jHT|) and solving it next with respect to Y: 

^ I Hi + 2ER-+y{i^ + 2) - i?-+V(z/ + 1) + /3/?2(.+2)/(^ + 2)2] 

^ - ^m^) 

It is evident from equation (|32p. that incoming and outgoing separatrices are symmetri- 
cal with respect to OR axis. Therefore we can restrict our analysis to one of them, e.g. 
to the upper separatrix Y+. First of all, let us note, that in the point 0) separa- 



trix y+ forms with OR axis a positive angle a = arctg^J{R2 — Ri)'^{Ri)/ (aRi^'^^. The 
above formula arises from the linear analysis of system ()26p in critical point Ai{Ri, 0). 
So Yj^{R) is increasing when R — Ri is small and positive. On the other hand, function 
G{R) = Hi + 2E i?^+2/(z/ + 2) - [2^)2 R'^+^/(u + 1) + f3 i?2(^+2)/(^y + 2)2] , standing inside 
the square root in equation tends to —00 as R ^ +00, because the coefficient at 

the highest order monomial i?2{!^+2) jg negative, while the index z/ + 2 is assumed to be 
positive. Therefore the function G{R) intersects the open set R > Ri of the OR axis at 
least once. Let us denote a point of the first intersection by A^lR^jO). The coordinate 
of intersection satisfies inequality R3 > R2. It can be easily seen by noting that G'{R) = 
-2R'' {PR''+y{u + 2)-ER + ^2) = -2R''P{R) and analyzing lim^^^- dY/dR. On the 
other hand, G'{R) < when R > R2, and this is suffice to show that lim^^^- dY/dR = —00. 
Due to the symmetry of the saddle separatrices, they intersect tangently in the point 
A^i^R^, 0), The result obtained can be formulated as follows. 

Theorem If u > —2 and > (3 Ri^^ , then system 1126]) possesses a one parameter 
family of periodic solutions and the homoclinic solution, formed by the tangent intersection 
of separatrices of the saddle point Ai 0). 

Thus we have shown, that system (j24|) possesses periodic and soliton-like invariant solu- 
tions. Let us note in conclusion that for some special case the integral standing at the RHS 
of the formula can be calculated explicitely: 



u = ±< Vo arcsm 



R+1 
2V2 



+ 72 log 



R-1 



3-R+J2O- 2(i?2 + 2i? + 3) 



-uJo\, (33) 
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where ujq = -\/2 vr — log2. This solution corresponds to the following values of the parameters: 
D = l = R^ = a = l, P = 1/2, z/ = 0, and E = 5/4. 

4 Conclusions. 

Thus we presented the effective anzatz resulting from the analysis of the methods put forward 
in OElHj and based on the following observation: the crucial element of any of anzatz-based 
method is the closeness of certain class of functions with respect to the algebraic operations 
and differentiating. This observation served as the main motive when we put forward the 
anzatz 0. Note, that it can be easily modified for obtaining the periodic solutions or, more 
generally, exact solutions that can be expressed as algebraic combinations of some special 
functions. 

The results of the last section show that the combination of the self-similar reduction and 
the qualitative analysis can deliver the exhaustive answer to the question on the existence 
of certain types of solutions within the given family. At the end of this section we also 
present the exact homoclinic solution, that was predicted earlier by the qualitative analysis 
but failed to be obtained within the anzatz-based methods. The reason of this is evidently 
linked with the fact that solution (jH^ is too comphcated to be expressed as an algebraic 
combination of the hyperbolic functions. Whether it can be expressed in terms of special 
functions is still an open problem, connected with the more general problem of finding out 
a fully universal ansatz. 
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